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Most large cities are spanned by more than one transportation system. These different modes of transport have 
usually been studied separately: it is however important to understand the impact on urban systems of the coupling 
between them and we report in this paper an empirical analysis of the coupling between the street network and 
the subway for the two large metropolitan areas of London and New York. We observe a similar behaviour for 
network quantities related to quickest paths suggesting the existence of generic mechanisms operating beyond the 
local peculiarities of the specific cities studied. An analysis of the betweenness centrality distribution shows that 
the introduction of underground networks operate as a decentralising force creating congestions in places located at 
the end of underground lines. Also, we find that increasing the speed of subways is not always beneficial and may 
lead to unwanted uneven spatial distributions of accessibility. In fact, for London - but not for New York - there 
is an optimal subway speed in terms of global congestion. These results show that it is crucial to consider the full, 
multimodal, multi-layer network aspects of transportation systems in order to understand the behaviour of cities and 
to avoid possible negative side-effects of urban planning decisions. 
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INTRODUCTION 

In the last decade, urban transportation networks have been largely studied by means of spatial networks [I] in 
order to understand aspects of urban systems and their evolution. These studies comprise the morphology of street 
networks [2-17], their evolution [18-21], and their relationships with socio-economical indicators [22, 23]. In parallel, 
there are also studies on the structure of subway networks [24-26], their evolution [27], and their robustness [28, 29]. 
However, these networks are not independent, and the important result in [30] showed that coupling between networks 
can be critical and can affect the global behaviour of a system. It is in this context that multilayer (or multiplex) 
networks [31-34] are studied and provide the convenient conceptual framework. A few recent studies considered the 
impact of the multilayer structure on various general processes [35, 36], and specifically in the case of transportation 
networks [37, 38]. Multilayer networks offer a good theoretical framework for understanding how interconnected 
transportation networks are shaping cities and how they may affect their operation. Moreover, given the increasing 
interest in urban systems, an empirical study of the impact of the multilayer structure of transport systems on mobility 
appears both crucial and timely. 

In this study we consider the mutually connected underground and street networks in the large metropolitan areas 
of Greater London and New York City, and explore how their coupling affects their global properties. In particular, 
we analyse the effect of varying the subway speed and show that increasing it can lead to unexpected counter-effects. 
Our analysis focuses on three main network features and findings: i) the behavior of quickest paths at the city scale; 
ii) the local outreach and the urban horizon; and iii) the spatial distribution of betweenness centrality. It is important 
to stress that studies of urban transportation networks have important implications for urban policies and private 
investment, and in general play an important role in the urban planning chain. In fact, inter-modality transportation 
efficiency and simulations have been extensively studied in the transportation engineering literature [39], where the 
typical supply-demand approach prevails but where the analysis of topological properties of networks is almost wholly 
neglected and where the different transportation modes are often treated separately. One goal in this study is to 
shift the focus onto this topological coupling aspect of transportation network design: we show this to be extremely 
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relevant, and suggests that the multilayer network view of these systems should be integrated into elaborated models 
of urban planning. 


DATA AND NETWORK CONSTRUCTION 

Using data from Open Street Map [40], we construct both the street and the subway networks for London (UK) and 
New York City (USA). We downloaded data on street networks and underground in geo-referenced vectorial format 
from Open Street Map, which contains detailed streets and rail tracks networks, including train depots and double 
tracks. (The rationale behind the geographical extent of these networks is to include the full underground systems 
and surrounding street networks.) In addition, a series of automatic and manual topological cleaning operations were 
needed in order to extract consistent and usable graphs.The size and geography of the two cities is clearly different 
as we can observe it in Fig. la,b. 





P 


FIG. I: (a,b) The spatial extent of the two metropolitan areas considered here. Note that the Greater London area 
is not covered by the underground system, in contrast with New York where most areas are connected by the 
subway, (c) Distribution of normalised quickest path times (computed in the multilayer system) zn{Tm) = 

y/Var(Trn) 

(d) The quantity (qms) averaged over all nodes as a function of (3 (the error bars indicate here the dispersion around 
the average). The average ratio of travel times with and without the subway layer is typically of order 0.5 and does 
not vary much with f3. (The figures Ia,b were created using vectorial map extracted in ArchMap environment and 

assembled with Adobe Illustrator.) 

We thus obtained the weighted graph Gg = (14, Ug) of the connected street network in its “primal” representation, 
with nodes being street junctions and edges representing the street segments connecting them, and the weights given 
by the street length W = streetlength{m). Similarly, we obtained the connected underground network Gu = (14,14) 
with nodes representing underground stations and links connecting successive stations on the same line, and weighted 
by the length of the line segment W = linklength[m). From a theoretical point of view, the interdependent or 
multilayer [31] network, Gmuiti is defined as the union of these two networks. Here we have subway stations and road 
intersections that we consider to be different nodes. Underground stations are accessible from more that one access 
on the street, but for the sake of simplicity we construct the multilayer network by connecting each underground 
station to its closest street junction only (a simplification that won’t change the structure of quickest paths). In order 
to create the adjacency lists we used a combination of Python scripts, ArchMap geo-tools (ArchGIS 10.2) and ad hoc 
manual corrections. Tools have been set to remove links redundancies, to correct the topology of the networks, and 
to create the proximity matrix between street nodes (street junctions) and underground nodes (stations). The scripts 
have been corroborated with a full check of the data and corrections of the topology in editing sessions under ArchMap 
environment (the computation of the various statistical measures have been done under Python environment using 
NetworkX library and the maps have been produced using ArchGIS 10.2). 
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THE GENERIC NATURE OF QUICKEST PATHS 

New York is composed of two large and almost disconnected components with the underground systems covering 
a similar spatial extent and carving-up the different boroughs. London instead presents - at a large scale - a typical 
radiocentric urban structure with the underground systems connecting satellite districts and peripheries to the urban 
core. Differences both in size and geography between these cities are also reflected by basic network descriptors shown 
in Table I. 


Case 

N 

M 

Total length (km) 

■jgeo , . 

hj (km) 

-rtopo 

Hj 

DiamP’^° (km) 

Diand°^° 

\/A(km) 

Node density (km ^) 

London streets 

324536 

427920 

34493.73 

25.83 

178.16 

89.31 

368 

73.68 

59.78 

London tube 

263 

296 

385.98 

18.55 

14.26 

60.3 

42 

- 

0.048 

London multi 

324799 

428479 

34886.52 

25.78 

96.16 

89.27 

288 

- 

- 

New York street 

68417 

112827 

12153.81 

17. 94 

106.64 

55.18 

278 

34.93 

56.07 

New York subway 

454 

489 

416.12 

18.87 

19.10 

57.28 

62 

- 

0.37 

New York multi 

68871 

113770 

12579.44 

17.91 

54.45 

55.18 

205 

- 

- 


TABLE I: The number of nodes N, the number of links M, the cost defined as the total length of all links, the 
average geographical and topological shortest path length, (km) and respectively; the maximum 
geographical and topological shortest path length, (km) and D*°p° respectively; the density of streets Ps = ^ 

and the density of underground stations Pu = ^, where A is the area of each city. 

For both cities, the (spatial) diameter of the multiplex is essentially dominated by the street network. We also 
observe that the topological diameter of the multiplex is lower than the street layer, thanks to the subway structure 
allowing for topologically shorter paths. The efficiency of the subway is however also due to its speed which is in 
general larger than that of overground modes such as private cars, taxis, or buses. In order to reflect this, we introduce 
a parameter 0 < /3 < 1 that describes the ratio of speeds in both systems, similarly to the theoretical analysis proposed 
in [37]. This parameter (3 measures the travel cost in time units associated to the underground links. This means 
that the number of time units taken to traverse an underground link of length I meters is (31, which is ^ times faster 
than the time taken to traverse the same length on the street network. Thus a smaller j3 corresponds to a faster 
underground speed, as compared with the speed on the street network. The introduction of this parameter allows us 
to study the properties of the multilayer system as a function of underground speed. Naively one could expect that 
the system as a whole will be more efficient for faster subways, but we will show here that it is not always the case 
and that in some cases we can observe an optimal value for j3. Finally (3 can be measured empirically, and we obtain 
for London ^London ~ 0.48 and a slightly larger value for NY (3nyc ~ 0.55. 

We denote by Ts{i,j) the travel cost (ie., the number of time units) of the quickest path between street nodes 
i,j G Us 5 and by Tm{i,j) the cost of the quickest path between i and j in the multilayer network (ie., a path which 
can traverse both street and underground links). The normalised quantity 


2/3 (r^) 



^Var{Tra) 


( 1 ) 


displays a behaviour that is roughly constant for (3 larger than 0.2 — 0.3, as shown in Fig. Ic, demonstrating that the 
effect of /3, is essentially contained in the average and variance of Tm- This is a rather surprising result, given that 
the two cities display many geographical and structural differences. The cost between nodes i and j can be 

written as 


eeP{i,j) 

where the sum is over all links e that belong to the quickest path P{i,j) and where T(e) is the cost on this link. (We 
neglect inter-modal change costs in this simple argument.) If the path is long enough, and if the random variables 
r(e) do not display long-range correlations and are not broadly distributed, the central limit theorem applies and the 
distribution of the Tm follows a Gaussian distribution in a certain range. There are obviously deviations observed 
for small values of (3 coming from the fact that the paths’ durations become very heterogeneous depending on the 
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proximity of their origin or destination to subway stations. In this respect, a very high relative subway velocity 
enhances spatial differences in the city and may lead to an uneven distribution of accessibility, a fact that will be 
confirmed below with the local outreach analysis. 

We also compute the average ratio between the travel costs from i to other street nodes through the multilayer 
network and through the street network, defined as 


Qms (f ) 


1 

AW 


E 




( 3 ) 


where Ns is the number of street nodes. The larger this ratio, the larger the effect of the underground on travel costs. 
We see in Fig. Id that typical values are of order 0.5 for both cities and that the effect of /3 is rather weak: a decrease 
from /3 = 1 to /3 = 0.5 leads to a decrease of {qms) of order 20%. In addition, it seems that the effect of subways in 
London is less important than in New York, which is probably due to the lesser extent of the subway in the Greater 
London area. 

A central quantity for describing the importance of inter-modality is given by 


j) 


—multi 




( 4 ) 


where Cij is the total number of shortest paths between i and j (using either one or two networks), and the 

number of paths using edges of both networks at least once. It characterizes the importance of multi-modality for the 
path from i to j. If we sum over all possible destination nodes j, we can quantify the added value of the interlayer 
coupling to the reachability of nodes, and obtain the interdependency [37] of a street node z € 14 defined as 


A(z) 


1 


E 


—multi 




( 5 ) 


(Note that a similar measure has been used in the transportation design literature under the name of inter-modal 
connectivity [39].) In order to understand the effect of scale on the interdependence, we also define the interdependence 
profile as 


Qx{d) 


1 

N{d) 


E 

de.{i,j)=d 


( 6 ) 


where de{i,j) is the Euclidean distance between i and j and N(d) is the number of pairs of nodes at Euclidean 
distance d. In Fig. 2a we show the average interdependence among all street nodes as a function of (3 and the 
resulting interdependence profile Fig. 2b. 

We see from these figures that, in both cities, the existence of the underground has a very large impact. For 
example, for (3 = 0.8 we obtain A around 0.7, meaning that even when the underground is only 1.25 times faster than 
the street network, already about 70% percent of the quickest paths are going through the underground. A slight 
decrease of /? for (3 close to one thus has a large impact on structure of quickest paths, while for smaller values of 
/3, improving the subway speed does not bring a significant improvement of the structure of quickest paths. In both 
cities there is a sharp increase in A for small Euclidian distance, meaning that already for relatively short trips, it is 
worth “hopping on” to the underground. (Note that we neglect here waiting, walking, and connecting times which 
can significant [38].) The slope of the interdependence profile at small deuci — 0 is increasing as j3 is decreasing, 
suggesting that a slight increase in the underground speed could make the networks highly interdependent even at 
very small scales. 

Both cities therefore display a remarkably similar behaviour over all these interdependency-related quantities (in 
particular, see Fig. 2(b)), suggesting here again a possible common behaviour for multiplex transportation network 
in cities. While further studies are needed in order to substantiate a claim of “universality”, our results point to the 
possible existence of some kind of a statistical law of large numbers that applies to quickest paths in multiplex urban 
transportation networks. 
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FIG. 2: (a) Average interdependence, (A) as a function of p. (b) Normalised interdependence profile computed for 
different values of /?. Both cities exhibit a similar behaviour despite very different geographical structures. 


We note that it is not trivial that the central limit theorem applies here and it doesn’t mean that the network 
topology is irrelevant. The fact that we can sum a large number of quantities, which are uncorrelated (a necessary 
condition for the central limit theorem to apply) comes from the specific structure of these transportation systems 
(spatial constraints for example certainly play an important role). In addition, more complex quantities (such as the 
interdependence for example) also display a large level of similarity for the two cities, a fact that cannot at this stage 
be simply related to a central limit theorem. These different results point to the potentially useful fact that actually 
few parameters seem to govern the behavior of these quantities, which could lead to many useful simplifications in 
more elaborated models that contain a large number of parameters. 


LOCAL OUTREACH AND THE URBAN SPATIAL HORIZON 

The presence of a transportation mode such as a subway affects the overall performance of a city in terms of 
efficiency of transport and the accessibility of certain locations, but also has an important impact on how pairs of 
locations are connected. In order to measure this effect, we define the spatial outreach of a street node i G Ustreet as 
the average Euclidean distance from i to all other street nodes that are reachable within a given travel cost, t: 

= N(t) ^ de{i,j) (7) 

where de{i,j) is the Euclidean distance between node i and j, and N{t) is the number of nodes reachable on the 
multilayer network within a given travel cost r. In Eig. 3 we show the average local outreach as a function of the travel 
cost threshold r, which displays a non-linear behaviour due to the different speeds achievable in the two transportation 
modes. This provides support for a general effect already known: for longer trips, faster transportation modes are 
used (see for example [38] for the UK case). 

For New York, the unit of time is given by the average car speed on the street network which is of order 15.6km/h 
(see for example [41]). Rescaling r by this velocity, we then obtain an effective maximum speed (for /? = 0.1) of 
order 30km/h (close to the 28km/h discussed in [42]). For London, the same calculation with an average car speed of 
16km/h (see [43]) yields an effective maximum speed of order 21km/h. (This difference in speeds is due to the areas 
considered, as New York is almost entirely covered by the underground network.) 

As shown is Fig. 4,a,b,d,e as P decreases, the nodes having a high local outreach are concentrated close to under¬ 
ground stations where the underground is the most accessible, and the graph consisting of high-outreach nodes (red 
nodes on the map) becomes less fragmented. 

In other words, as the underground becomes faster, a continuous area of high-outreach nodes emerges (the com- 
mutable zone) in the city centre and around the nodes of the underground network, implying that a person can travel 
from this area to faraway places (large Euclidean distance) at a small travel cost r. The location of this highly 
accessible zone cluster from a dispersed configuration, as in Eig. 4d, to a centralised one, as in Eig. 4a, which shows 
a centralisation effect due to the accessibility provided by the underground. The dispersion of the local outreach also 
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a 


London 


b 


New York 




FIG. 3: Average 


local outreach {L) normalised by the square root of area, for London (a) and New York (b). 




London 



P 

New York 



P 


FIG. 4: Maps showing for London the spatial distribution of the local outreach defined in Eq. 7 for a travel cost 
r = 2000 and speed ratio ( a) /?=0.1, (b) /3=0.7, (c) /3=0.4 and (d) /3=1. As the underground’s speed is increasing 
compared to the speed of the street network (ie., decreasing P) the nodes having a high local outreach are 
concentrated at the center along the underground network, where the underground is most accessible. Theses figures 
were created using ESRI ArchMap 10.1 and Adobe Illustrator, (e, f) The graphs show the Gini coefficient of the 

local outreach versus /3 for both cities. 


displays a very interesting result demonstrated by its Gini coefficient Gl € [0,1]. Indeed, in Fig. 4e,f we see that for 
both cities for /3 > 0.5 the accessibility is distributed almost uniformly amongst all the places in the cities, while for 
smaller /3 (faster underground) the shift to an uneven distribution of accessibility is clear. This result suggests that 
transportation policies that focus on increasing the speed on a single travel modality may lead to the undesirable 
spatial heterogeneity in the accessibility of different locations. 

We show in Fig. 5 the probability that the outreach is larger than a certain fraction aL of the size of the city, and 
we observe the existence of a threshold ac- 

The existence of a threshold less than one means that, for given values of /? and r, there is a maximal value Lm for 
the outreach. We can estimate the value of Lm by using a simple argument: the maximum value is reached when the 
path is “essentially” made on the quickest transportation mode, the subway. This transportation mode has a velocity 



















7 


0.01 


0.1 

a 





0.01 


0.1 



20000 

T/P 


FIG. 5: (a)-(b) Fraction of nodes with a local outreach larger than aVA for London (a) and New York (b). (c)-(d) 
when both axis are rescaled according to Eq. 8, we observed a data collapse with all curves collapsing on a single one. 


given by and the probability that a station is within reach (in circle of radius do corresponding to the typical 
walking distance to reach the subway) is 


p = pu-n-dl (8) 

where = N^/A is the density of subway station {A = L? is the area of the city and is the number of subway 
stations). The maximal outreach Lm is then given by 


L2 


- • u ,2 ^ 

= WwTranw'! 


( 9 ) 


and Oc = LmiL is then given by 


(y.c 


Nu 

L3 


Tid^v 


P 


( 10 ) 


This last equation shows in particular that the quantity a^.L'^/Ny_ should increase linearly with t//3, with a constant of 
proportionality depending on geometry the city, and we observe that this scaling in is agreement with simulations (see 
figures 5c, d). In particular, we see that the slopes for London and New York are different: the ratio of the constant 
pre-factors is about 10, suggesting that the subway system in London is more efficient in terms of the outreach that 
can used as a measure of the “urban horizon”. 


THE GEOGRAPHY AND DISTRIBUTION OF URBAN CENTRALITY 

Betweenness centrality (BC) [44] is one of the important quantities in complex networks, and in street networks 
in particular [8]. It quantifies the importance of a node as being the amount of traffic going through it, assuming 
uniform demand where the traffic between all pairs of nodes is the same. This quantity is very relevant in urban 












systems: in particular, it is correlated with the locations of shops and other micro-economic activity [22, 23], urban 
growth [19, 20], and land-use intensity [45]. 

In the case of car traffic and congestion, the absence of detailed traffic models or mobility data leads us to use the 
BC in order to identify the potentially congested locations and the effects of spatial structure on the shortest path 
structure. Even if we know that the assumptions used in the BC calculation can lead to some inacurracies [46], it 
is the simplest proxy that contains some level of information about real traffic. We thus explore in this section the 
spatial distribution of BC in the street network and how it is affected by the underground system. The BC of a street 
node V € 14 in the street network is defined as 


bcs{v) 


1 

(iV,-l)(iV,-2) 


E 


^street (^) 
^street 


( 11 ) 


where is the number of quickest paths between i and j in the street network, of which {v) goes through 

street node v. Similarly, we define the betweenness centrality of a street node G 14 in the multilayer network as 


bcyn{v) 


1 

(TVs - l)(iV, - 2) 


E 




( 12 ) 


where is the number of quickest paths between i and j in the multilayer network, of them goes through 

street node v. 

We can then observe how the parameter (3 impacts the mobility distribution and the geography of potentially 
congested areas. The maps in Fig. 6(a,b,c,d) show the BC spatial distribution for both cities computed on streets for 
/3 = 1 (a,b) and (3 = 0.1 (b,c). 

These maps clearly display a dramatic change in the spatial distribution of central places when introducing an 
underground system, shifting congestion from internal street routes and bridges to inter-modal places located at the 
terminal points of the underground networks, which presumably are used as entry/exit gates for suburban flows to 
reach core urban areas. Remarkably enough, in both cities, these places are located in urban areas that do not 
overlap with the underground system, thus possibly creating congestion in unexpected places. In other words, the 
introduction of underground networks operate as a decentralising force creating congestion in places located at the 
ends of underground lines and not, for example, in the city centre as one might expect referring to classical results 
on rewiring processes for chain or lattice networks [1] in which BC is correlated to the distance to the gravitational 
centre. The statistical dispersion of BC can be measured by its Gini coefficient and also suggests that congested places 
always become more critical in the system as (3 decreases. In fact, as shown in hgure 7, the Gini coefficient of BC 
increases as the underground becomes more efficient (faster, decreasing /3), meaning that a larger fraction of quickest 
paths use it; and the BC distribution is less homogeneous, making the system more fragmented and less resilient. 

Examining BC Gini as a function of (3 and the interdependency A in London (Fig. 7), we observe a non-trivial 
optimal value for /3 for which flows are the most homogeneously distributed across street junctions. This minimum 
value for the Gini also corresponds to a minimum of the total congestion (see the appendix Total congestion). In New 
York, however, there seem to be room for small (3 and small congestion and the absence of a non-trivial optimum for 
New York suggests (as discussed theoretically in [37]) that - surprisingly - it has a more marked monocentric aspect 
than London. In other words, the congestion in central places in New York is so large that introducing an efficient 
subway system is always better, even if it creates congestion at other points. Remarkably, these results on BC and 
on the existence of an optimal point are thus in agreement with a recent theoretical model of coupled transportation 
networks, where - depending on the distribution of trip targets - two regimes were observed: one in which the optimal 
coupling is trivially the maximum, an another one where a non-trivial optimal coupling exists [37]. 


DISCUSSION 

We have considered the effect of the coupling between two transportation layers on various quantities and we 
can summarise our results as follows. For quantities relating to quickest paths (interdependency, average quickest 
path duration), we observe a remarkable similarity between the two cities considered, suggesting the possibility of a 
universal behaviour requiring further study. This universality might originate in the fact that the quickest path can be 
seen as a sum of random variables, which inevitably leads to some sort of central limit theorem. This seems to be the 
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FIG. 6: The spatial distribution of BC on the New York (top) and London (bottom) street network for different 
values of p. We observe a clear crossover from congested road locations for /? = 1 to “focal points” of the 
underground system for small /3 (The maps have been obtained using the classical interpolation method Inverse 
Distance Weigh (IDW) on the street junctions with Z = BC). This figure was created using ESRI ArchMap 10.1 and 

Adobe Illustrator. 


case for the probability distribution of the quickest path time duration, that (once normalised) is a universal function 
for a reasonable range of subway speeds. More involved quantities such as the local outreach and the urban horizon 
also display a simple common behaviour that cannot be recovered using a back-of-the-envelope argument about the 
quickest path. This possible universality suggests that few parameters seem to govern the behavior of quickest paths, 
which could lead to many useful simplifications in more elaborated models that contain a large number of parameters. 
More data on more cities is however needed in order to validate this universal idea and to understand its origin. 

We also studied the impact of the coupling on the spatial distribution of the betweenness centrality. We observe 
results in agreement with previous theoretical findings, with in particular the existence for London of an optimal sub¬ 
way velocity in terms of congestion. These results on spatial distribution of centralities can be also understood in the 
framework of another study showing how decentralising housing in London leaves room for commercial activities [47]. 
Even if the direct causal link between land-use change and transportation network evolution is still not clear [48], our 
results seem to go in that direction. It would be interesting to confront our results with realistic models used by the 
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FIG. 7: We represent here for London (a) and New York (b) /3 


and interdependency A as a function of the dispersion 


of the multilayer BC (measured by the Gini coefficient Gb for the BC). 


transportation community, in particular when the subway speed is modified. More generally, we believe that more 
empirical studies are however needed in order better to understand the complex coupling between land-use and the 
structure of multimodal networks. 

It thus seems to be clear that it is important to consider full multimodal, multi-layer network aspects in order 
to understand the behaviour of an urban transport system - and thus understand the effects of transport on other 
features of interest. Even if these studies are still very theoretical, they show convincingly that reasoning with only 
one transportation mode can be extremely misleading, and that policy-makers cannot limit themselves to a single 
aspect of an urban system without risking making decisions that are locally correct but globally wrong. 


DATA ACCESSIBILITY 

Original network data is available on Open Street Map [40]. The cleaned and topologically corrected networks are 
also freely available [49]. They include ArchMap shape file (.shp) containing streets and underground networks with 
their adjacency lists. These files can be opened on any GIS platform. 
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APPENDIX: TOTAL CONGESTION 


We can choose as a measure of congestion the total (or equivalently the average) time spent on roads. We use the 
standard Bureau of Public Roads function [50] that gives the time to go from i to j (separated by a distance d{i,j) 
on a road with traffic T(i,j) and capacity c 


d{i,j) 


1 + 




(13) 


where v is the average velocity on this road and where the exponent v is usually large « 4 — 5. In our model, the 
proxy for the traffic is given by the BC and the total congestion Ct is thus approximately equal to (up to a constant 
factor independent from j3) 


Ct oc ^-^{{bcr) (14) 

V 

(we assumed that the length distribution of segments is peaked around its average (d)). In the figure 8, we plot this 
total congestion - normalized by its value at /3 = 1 - as a function of /?. We observe on this hgure a clear optimum 



FIG. 8: Total (normalized) congestion computed for v = 4 (left) and maximum congestion (right) as a function of (3. 


which corresponds to the minimum of the Gini coefficient. Congestion is actually very sensitive to large values of 
the traffic and since the maximum has also a minimum (Fig. 8) we can indeed expect that the larger the traffic 
heterogeneity and the larger the total congestion. 
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